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Differential  multiscale  modeling  of  chemically  complex  materials  under  heavy  deformation:  Biological,  bioinspired 
and  synthetic  hierarchical  materials 

ABSTRACT 

This  research  was  focused  on  modeling  and  design  of  high  stress  and  impact  mitigating  structures,  utilizing  nanoscale  patterning  and 
hierarchical  biomimetic  concepts.  The  eventual  goal  is  to  create  heterogeneous,  hierarchical  designs  for  thin  coatings  and  bulk  materials, 
capable  of  providing  enhanced  ability  to  mitigate  high  rate  impact  and  deformation.  The  potential  of  a  structure  to  mitigate  impact,  large 
stress  and  large  deformation  is  characterized  by  (i)  the  ability  of  the  material  to  dissipate  energy  under  high  rate  deformation,  (ii)  the 
resistance  to  brittle  fracture  by  crack  formation 

under  high  rates,  and  (iii)  the  ability  to  redistribute  load  underneath  a  thin  external  coating  film.  To  achieve  this,  our  efforts  are  centered  on 
the  development  of  a  holistic  atomistic  based  core  model  of  the  deformation  and  fracture  mechanisms  of  thin  nanostructured  coatings.  Using 
atomistic  simulation,  we  study  the  behavior  of  nanostructured  composites  under  heavy  impact  loading,  incorporating  different  material 
combinations  that  are  coupled  in  various  arrangements,  at  different  length  scales,  arranged  in  a  hierarchical  pattern.  The  material 
combinations  feature  divergent  characteristics  such  as  hard-soft  or  brittle-ductile,  since  mixture  of  materials  with  disparate  properties  are 
often  found  in  Nature’s  toughest  and  mechanically  most  robust  materials,  used  to  provide  protective  surfaces  (e.g.  in  seashells,  bone,  spider 
silk).  We  demonstrated  the  development  and  application  of  such  material  design  paradigms  in  studies  of  bone,  silk,  collagen  (and  similar 
materials),  enabled  through  the  development  of  multiscale  models.  Our  work  has  pushed  the  frontier  of  biomechanics  and  biomaterials 
modeling  to  enable  a  bottom-up  perspective  of  key  issues  that  define  robustness,  strength  and  adaptability  of  biological  and  biologically 
inspired  mechanically  relevant  materials  for  numerous  applications. 
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0  Differential  multi-scale  modeling  of 

chemically  complex  materials  under  heavy 

it .  deformation 

Motivation 
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-  Need  for  lightweight  shock-mitigating  coatings 

-  Develop  inspired  from  biological  structural 
materials  (bone,  nacre),  combine  disparate 
material  properties 

-  Challenge:  To  optimize  hierarchical 
nanostructured  material  using  a  bottom-up 
atomistic  modeling  approach 


NanDstnjcturad  coating  Otnm} 


•  More  details  about  the  scientific 
challenge 

-  Limited  understanding  of  hierarchical 
nanostructures  and  of  integration  of  structure 
and  material  in  biological  materials,  no 
theoretical  constitutive  models 

-  Prevents  from  developing  engineering 
materials  with  similar  hierarchical  features 


Macroscale  structural  engineering 


This  research  project  is  focused  on  the  analysis,  modeling  and  design  of  impact 
mitigating  coatings  and  thin  films,  utilizing  nanoscale  patterning  and  hierarchical 
biomimetic  concepts.  The  goal  is  to  create  heterogeneous,  hierarchical  design 
suggestions  for  thin  coatings  capable  of  providing  enhanced  ability  to  mitigate  high 
rate  impact  and  deformation.  The  key  new  element:  Use  of  atomistic  based 
multi-scale  modeling  to  explore  the  opportunities  and  potential  of  including 
the  atomistic,  molecular  and  mesoscale  structure  in  the  engineering  design 
space  by  an  integrated  treatment  of  structure  and  material  (M.J.  Buehler  and  Y. 
Yung,  Nature  Materials,  2009;  M.J.  Buehler  and  T.  Ackbarow,  Materials  Today, 

2007,  Ackbarow  et  al„  PNAS,  2007). 

The  potential  of  a  structure  to  mitigate  impact  is  characterized  by  (i)  the  ability  of  the 
material  to  dissipate  energy  under  high  rate  deformation,  (ii)  the  resistance  to  brittle 
fracture  by  crack  formation  under  high  rates,  and  (iii)  the  ability  to  redistribute  load 
underneath  a  thin  external  coating  film. 

To  investigate  these  properties,  our  efforts  are  centered  on  the  development  of  a 
holistic  atomistic  based  core  model  of  the  deformation  and  fracture  mechanisms  of 
thin  nanostructured  coatings.  Using  atomistic  simulation,  we  study  the  behavior  of 
nanostructured  composites  under  heavy  impact  loading,  incorporating  different 
material  combinations  that  are  coupled  in  various  arrangements,  at  different  length 
scales,  arranged  in  a  hierarchical  pattern.  This  insight  is  then  used  to  develop 
design  suggestions  for  new  nanocomposite  structures. 


2 


$  Biological 


Biological  systems  universally 
feature  active  components; 
e.g.;  bone  remodels  according 
to  changes  in  loading,  cells 
change  their  shape  under 
different  environments  or 
stress 

Mechanisms:  Sensing  of 
environmental  cues  followed 
by  gene  regulation/activation, 
which  changes  the 
hierarchical  structure  of 
materials  ("mutability") 


materials  are  adaptable, 
mutable... 


Process 
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New  blood  vessels  form  under  cyclic  strair 
to  repair  wounds,  provide  nutrients  where 
needed 
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Mechanotransduction  mechanisms  (e.g,  angiogenesis} 

Bushier  and  Yung,  Nature  Materials ,  2009  y  Yung.  M.  Buehler,  D,  Mooney,  ef  at.,  PNAS,  2009 

Challenge:  Time-scales  on  the  order  of  minutes,  hours,  days  3 


M.J.  Buehter  MfT 


0  Differential  multi-scale  modeling  of 
chemically  complex  materials  under  heavy 
deformation 

•  Dbjective 

-  Develop  nanostructured  material  that  resists 
fracture  (cracking)  and  redistributes  impact 
energy  into  plastic  dissipation  over  large  volumes 

•  Approach 

-  Atomistic  based,  bottom-up  multi-scale  simulation 
combined  with  theoretical  modeling,  to  quantify 

the  role  of  length-scales  and  material 
hierarchies 

•  Uniqueness  of  approach 

-  Our  approach  extends  current  state  of  the  art 
(phenomenological  models/approaches)  by 
providing  a  fundamental  description  of  material 
behavior  under  impact  (cracking,  plasticity,  shear) 

-  Provides  opportunity  to  design  material  from 

_ bottom  up _ 


Hierarchical  structures 
with  characteristic 
features  at  multiple 
scales 


This  slide  summarizes  the  objective,  approach  and  uniqueness  of  the  approach. 
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Differential  multi-scale  modeling  of 
chemically  complex  materials  under  heavy 

deformation  m.j.  Buehierwu 


Time 

scale 


Slide  shows  overview  over  multi-scale  simulation  paradigm.  The  integration  of 
several  numericals  that  operate  at  distinct  length-  and  time-scales  enables  us  to 
traverse  from  fundamental  chemical  scales  (atom-atom  interactions)  to  larger  scales 
that  are  relevant  for  applications  (that  is,  for  measurement  of  material  performance). 
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Why  mimic  biological  structures? 

M.J.  Buehler,  MIT 
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Banana-cuwe  of  stiffness  and  fracture 
toughness  of  materials1-  How  do  we 
mimic  behavior  on  the  right  figure? 
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Stiffness  and  toughness  of 
biological  materials2  3 


1000 


1.  Ashby.  Materials  Selection  in  Mechanical  Design  (1999);  2,  Fratzl  et  al.  J,  Mater  Chem,  14.  2115 
(2004);  3.  Ashby  et  al,  Proc.  R.  Soc  Lond.  A  450,  123  (1995) 


The  structure  of  biological  materials  such  as  nacre,  bone  and  their  hierarchical 
organization  leads  to: 

a)  not  following  the  usual  banana-curve  for  synthetic  materials  of  strength  vs 
toughness  when  we  compare  materials  with  maximal  toughness  and  maximal 
strength, 

b)  not  following  rule-of-mixtures  for  composite  materials,  and 

c)  incorporating  best  properties  of  constituent  materials. 


For  engineering  materials,  shown  in  the  left  figure: 

1  -  Conventional  engineering  materials  (e.g.  ductile  metals)  are  tough  but  lack 
strength 

2  -  Conventional  strong  engineering  materials  (e.g.  ceramics)  show  high  strength 
but  lack  toughness 

3  -  Hierarchical  nanostructured  materials  enable  to  reach  high  strength  and 
high  toughness  (illustrates  the  potential  of  utilizing  nanostructures  and  hierarchies) 
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Mimicking  biological  structures 


M.J  Buehler,  MIT 

One  level  design  (nanostructural  design):  universality  only  - 
maximize  utilization  of  sub-optimal  building  materials:  not  able  to 
have  different  properties 

Hierarchical  design:  universality  PLUS  diversity  for  multiple 
properties  (address  additional  design  needs)  — 


Strength/ 
"►  toughness 


1  Hierarchy 


2  Hierarchies 


Summary  of  basic  approach  in  designing  materials  with  multiple  hierarchy  to 
optimize  strength  and  toughness. 
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Hierarchical  protein  structures 


Effect  of  hierarchies 

-  Enable  combination  of 
disparate  properties 
(strength  and  robustness) 

-  Overcome  limitation  of 
many  engineering 
materials 


M  J  Bushier,  MIT 

—  Alpha  helix 
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Same  number  of  AHs  available 
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Analyze  distribution  of  strength  and 
robustness 

Arrange  differently  in  hierarchies 

Ackbarow  and  Buehter,  Nanotechnology,  2009 


This  slides  summarizes  results  obtained  by  our  Hierarchical  Bell  Model,  showing 
that  hierarchies  are  a  material  design  strategy  to  combine  disparate  properties 
(here  strength  and  robustness)  and  thereby  overcome  the  limitations  of  many 
engineering  materials. 


This  first  study  combined  only  8  elements,  which  we  arranged  in  all  possible 
hierarchical  combinations  (published  in  Nanotechnology,  2009).  This  result  shows 
that  by  solely  changing  the  arrangement  of  alpha-helical  proteins,  one  can  control 
the  strength-robustness  performance  of  materials. 


8 


We  then  extended  the  study  to  include  32  elements,  which  we  arranged  in 
all  possible  hierarchical  combinations  (published  in  IJAM,  2009;  inaugural 
issue  and  cover  paper).  The  results  confirmed  our  earlier  findings,  and 
shows  an  even  more  pronounced  “inverse  behavior”  -  suggesting  that  by 
arranging  the  alpha-helical  building  blocks  into  hierarchies,  a  significant 
performance  increase  can  be  achieved. 

The  “Banana  curve”  (seen  in  conventional  synthetic  materials)  is  indicated 
the  plot  as  well. 


For  more  than  10,000  elements  (also  reported  in  the  I  JAM  paper),  a  very 
interesting  result  is  obtained.  First,  the  banana  curve  behavior  is  recovered. 
Second,  a  subset  of  all  geometries  shows  the  inverse  banana  curve 
behavior.  Interestingly,  only  2%  of  all  structures  give  the  inverse  behavior. 
This  means  that  most  random  arrangements  of  structures  would  lead  to  the 
banana  curve  behavior. 


This  result  is  important  since  it  suggests  that  it  is  crucial  to  control  the 
specific  nanoscale  arrangements  of  building  blocks  in  order  to  reach 
high  material  performance. 
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Summary  of  the  interplay  of  size  effects  and  hierarchical  structures.  The 

combination  of  size  effects  and  hierarchies  is  crucial  to  achieve  macroscopic 
superior  performance  to  combine  strength  and  robustness. 
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/s  structural  optimization  at  multiple  length  scales  (= hierarchies ) 
the  key  to  combine  disparate  material  properties? 


This  slide  illustrates  how  biological  materials  overcome  the  limitation  of  the 
“banana  curve”.  A  possible  explanation  may  be  the  existence  of  material 
hierarchies.  Our  model  (previous  slide)  illustrates  this  effect  at  a  fundamental, 
atomistic  and  molecular  level. 
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Transfer  to  metallic  nanocomposite 
analysis  &  design 


M.J.  Buehler,  MIT 


*  Transfer  to  metallic  nancomposite  analysis  &  design 

Bone:  Composite  of  rigid/plastic  materials 


Transfer  biological 
structures  to 
engineering  design 


jr-rtin 

* 


r  =  1 1 1  o] 


(b) 


Slides  shows  the  next  step,  the  transfer  of  the  results  obtained  from  studying 
biological  protein  materials  towards  the  design  of  metallic  nanocomposites.  The 
lower  part  of  the  slide  displays  our  model  system  consisting  of  a  ductile  matrix  and 
stiff  inclusion,  arranged  here  in  a  bone-like  geometry. 


Multi-scale  approach 


M  ,J,  Buehter,  MIT 


100  pm 


100  nm 


Atomistic  model  of 
nanocomposite 

Pull  atomistic  simulation 


(e  g.  EAM) 


1  pm 


—►  a y(w, — — ►  v!  v!  v  » 

Extract  constitutive  law  ■*■*■*■*■*• 

(various  design  parameters,  r  ■*■*■■*■  * 
platelet  width ,  etc.)  Mesoscale  model 

Parameter  studies  (e  g  (triangular  lattice  of 

elasticity,  flow  stress)  particles) 

Interparticle  potential  fitted 
io  constitutive  law 


Design  level 

(e.g.  gradient,  hierarchical  structures, 
topology);  structural  details  at 
micrometer  scale  and  beyond 
Optimization  algorithm 


The  development  of  design  tools  for  materials  with  multiple  levels  of  hierarchical 
assembly.  From  the  fully  atomistic  study  of  mechanics  of  bone-inspired 
nanocomposites,  we  obtained  strength  of  the  material  as  a  function  of  shape  and 
size  and  arrangement  of  the  composite  phases.  The  next  step  is  to  use  these 
nanocomposites  as  building  blocks  for  a  second  level  of  hierarchy.  The  aim  is  to 
create  multi-hierarchy  composites  that  can  improve  properties  such  as  strength  and 
toughness  simultaneously  through  the  use  of  design  at  several  length  scales.  The 
figure  presents  an  overall  multi-scale  approach  with  two  levels  of  discretization  to 
solve  this  problem.  A  coarse-grained  method  is  used  for  the  second  hierarchy  level. 
The  integrated  use  of  atomistic  and  coarse-grained  methods  will  enable  us  to 
explore  the  broad  range  of  structural  scales  relevant  to  initiate  a  new  paradigm  in 
multi-scale  structural  engineering. 
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Transfer  to  metallic  nanocomposite 
analysis  &  design 


*4i 


M.J.  Buehler,  MIT 


Systematic  design  analysis  of  bio-inspired 
nanocomposites 
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We  begin  our  investigation  with  the  analysis  of  a  single  hierarchy  material  -  representing  a  combination  of  a  soft  and  stiff 
phase,  forming  a  hetero-nanostructured  composite.  This  is  the  first  step  to  transfer  the  insight  from  biological  materials 
to  the  design  of  bioinspired  composites, 


Goal  of  this  study:  Can  we  tune  the  large-deformation  mechanical  properties,  and  if  yes,  how? 


The  first  important  contribution  of  this  work  is  the  development  of  a  model  system  that  enables  us  to  elucidate  the  effect  of 
nanostructural  arrangements  (geometry,  size,  spacing,  ..)  on  the  large-deformation  and  impact  loading  properties.  Based  on 
full-atomistic  simulations  and  a  multi-scale  integration,  we  have  then  identified  a  set  of  important  design  parameters, 
including  (i)  design  at  the  nano-level  (nanocomposite),  (ii)  arrangement  of  design  elements  in  hierarchies,  (iii)  interfacial 
properties  and  (iv)  environmental  effects  on  interfacial  properties.  These  parameters  represent  critical  factors  for  biomimetic 
design  of  impact  resistant  structures. 

For  example,  the  atomistic  simulations  reveal  which  shape,  size,  spacing  of  second  phase  platelets  produces  the  maximum 
strength,  and  what  changes  in  deformation  mechanisms  are  seen  across  these  geometric  changes.  A  representative  set  of 
results  will  be  reviewed  in  the  next  few  slides. 

The  main  finding  of  this  large-scale  atomistic  simulation  study  is  that  there  exists  an  optimal  size  of  the  platelet  inclusions  that 
leads  to  maximum  flow  stress  (N.  Broedling  et  al. ,  J.  Mech.  Phys.  Solids,  2007;  D.  Sen,  Acta  Mat.,  submitted).  The  result  is 
significant  in  two  ways:  (i),  it  shows  that  we  can  indeed  tune  the  large-deformation  mechanical  properties  and  optimize 
dissipation  by  tailoring  the  nanostructure,  and  (ii),  it  shows  that  there  exists  an  optimal  length  scale  at  which  the  flow  stress 
and  thus  dissipation  reach  a  maximum.  For  applications,  this  result  tells  us  the  optimal  particle  size  that  should  be  dispersed 
in  a  soft  matrix  materials  (for  the  example  of  Ni-AI  the  optimal  size  is  approximately  40  nm). 

We  have  further  discovered  the  critical  geometric  parameters  for  maximum  yield  stress  at  large  deformation  (D.  Sen  and  M.J. 
Buehler,  Acta  Mat.,  submitted).  The  studies  revealed  mechanistic  insight  into  the  nanoscopic  stress  distribution  and  its 
relation  to  the  toughening  mechanism.  This  result  is  significant  since  it  reveals  that  not  only  the  size,  but  also  the  details  of 
the  platelet  arrangement  is  significant. 

For  little  to  no  platelet  offset,  there  is  no  significant  strengthening  effect.  For  platelet  offset  larger  than  0.25  the  maximum 
effect  reaches  a  plateau  value. 

These  results  will  guide  laboratory  processing  techniques. 
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Transfer  to  metallic  nanocomposite 
analysis  &  design 


* 


M.J.  Buehler,  MIT 


Current  study:  Implemented  coarse-grained  model  to  reach  much 
larger  length-scales 

Constitutive  behavior  derived  from  atomistic  simulation,  enables  us  to 
reach  much  larger  scales 


We  have  developed  and  tested  a  mesoscale  coarse-grained  model  (similar  to  a 
finite  difference  method),  which  provides  a  bridge  from  atomistic  to  continuum 
length-scales.  We  are  currently  testing  the  performance  of  a  variety  of 
microstructural  arrangements,  where  each  phase  shown  here  consists  itself  of  a 
specific  nanostructure  (see  previous  slides).  This  provides  a  seamless  integration 
from  atomistic  to  continuum  and  the  tool  could  be  very  useful  for  design 
applications. 


16 


A 


4^ 


Multi-scale  coupling:  Coarse-grained  model 


kxE 


The  details  of  the  constitutive  law  for  the  coarse-grained  system  derived  from  fully- 
atomistic  simulations.  This  force  law  is  used  for  the  coupling  between  scales  of  the 
two  hierarchy  levels,  (a)  shows  the  typical  bulk  uniaxial  stress-strain  behavior  for  a 
metallic  nanocomposite  structure,  from  which  elastic  modulus  and  flow  stress  can 
be  derived;  (b)  shows  a  model  \  being  fit  to  this  loading  behaviour,  the 
proportionality  constants  (in  the  inset  equations)  depend  on  lattice  structure  for  the 
particle  model,  and  are  derived  using  the  Cauchy-Born  rule;  (c)  shows  a  potential 
application  of  the  mesoscale  model  in  a  2-phase  2-hierarchy  material  with  a  surface 
crack  and  material  microstructural  variation  around  it. 
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The  particular  metallic  nanocomposites  chosen  for  the  2nd  hierarchy  level 
composite.  We  use  a  soft  and  hard  phase  at  2nd  H  level  with  same  stiffness  but  a  10 
times  difference  of  yield  strength.  This  can  be  achieved  by  staggering  platelets  at 
the  1st  hierarchy  (nanocomposite)  level. 
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Inspiration  from  sea  sponge  skeleton  structur 


protein 


silica 


Structure  in  sea  sponge  composed 
of  silica-protein  composite 


A 


Objectives: 

a)  Minimize  depth  of  plastic  yield  in  interior  {for  a  surface 
crack) 

b)  maximize  plastic  zone  size  (toughness) 

c)  maximize  overall  stress  carried  (strength) 


Mesoscale  model  of  confinuous-laminated 
composite  desired  plastic 


zone  shape 


The  first  mesoscale  structure  we  study  is  a  laminated  composite.  Laminated 
composites  are  found  in  protective  skeletal  structures  in  several  biological  systems 
such  as  sea  sponges,  corals.  The  laminated  design  in  these  materials  is  critical  for 
imparting  damage  tolerance  to  the  overall  structure  (Seshadri  et  al.,  2002;  Chai  and 
Lawn,  2002).  Our  aim  is  similar,  to  maximize  fracture  toughness  to  surface  cracks 
while  minimize  their  plastic  deformation  in  the  interior  of  the  materia. 
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Preliminary  results  of  the  laminated  2-hierarchy  composite.  A  surface  crack  is 
created  under  remote  mode  I  loading.  The  volume  fraction  of  the  softer  phase  is  0.2 
and  constant  in  all  systems,  but  the  spacing  ratio  between  layers  is  varied.  The  red 
region  indicates  all  material  that  has  yielded  at  a  particular  loading.  The  effect  of 
reducing  the  spacing  between  soft  layers  lead  to  a  more  spread-out  plastic  zones, 
with  less  penetration  of  the  plastic  zone  in  the  interior  of  the  material. 
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AV  Material  choices  for  bone-inspired 

nanocomposite  structures  m.j.  Bushier,  mit 

•Soft  metal  (Al,  Cu)  matrix  -  hard  metal  platelets  (W) 

•Soft  metal  (Al,  Cu)  matrix  -  ceramic  platelets  (Al203) 

•Polymer  matrix  (PM M A/polystyrene)  -  CNT/graphene  platelets 


Polystyrene-graphene  nanocomposites  of  different  volume  fractions1 
1 ,  Stankovich  et  al,  Nature  442,  282  (2006), 


M  X  Buehterr  MIT 


New  concept:  Graphene  mechanics  in 
extreme  conditions 
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Characteristics  of  the 
polymer-graphene  interface 


M  J  Buehier,  MIT 


Morphology  of  the  graphene  sheets  in  PMMA/graphene  composite1 


How  strong  should  the 
adhesion  be? 


Fracture  at  the  polymer-graphene  interface1 


1-  Ramanathan  et  al,  Nature  Nanotechnology  3,  327  (2008), 


Persistent  wrinkled  nature  of  the  functionalized  graphene  sheets  within  the 
composite  provides  for  better  interaction  with  the  host  polymer  matrix.  Lower  figure 
shows  fracture  surface  topography. 
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Tearing  of  graphene  from  substrates 


M.J  Buehler,  MIT 

Tearing  mode  of  failure  at  graphene-polymer  interface  modeled  through  simpler 
system  of  graphene  tearing  from  model  substrate 


SEM  images  of  fragments  broken 
from  single  graphene  sheet 


Simulation  setup 


The  strength  of  graphene-polymer  interfaces  is  crucial  to  the  overall  load-carrying 
and  failure  mechanisms  of  the  composite.  A  the  interface  between  PMMAand 
graphene,  a  significant  failure  mode  may  be  tearing  of  the  graphene  sheets  from  the 
interface.  Here  we  look  at  the  influence  of  the  adhesion  strength  on  the  tearing 
mechanics.  Figure  on  the  left  shows  fragments  of  graphene  torn  from  larger 
graphene  sheets,  showing  a  tapering  behavior.  On  the  right  is  the  design  of  the 
computational  experiment  to  study  this  behavior  using  molecular  dynamics. 
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The  torn  sheet  is  seen  to  taper  off  at  all  adhesion  strengths.  We  have  measured  the 
angle  of  taper  as  a  function  of  strength  of  adhesion,  (a)  shows  plot  of  sine  of  the 
angle  of  tear  vs  square  root  of  the  adhesion  energy;  (b)  shows  predicted  angle  of 
tear  for  graphene  adhesion  on  graphite  as  a  function  of  number  of  sheets  torn  off 
together. 
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Shows  the  atomic  level  virial  stresses  for  the  graphene  sheet  (a)  and  (b)  show  the 
out-of-pane  shear  stress  siggZX  for  weak  and  strong  adhesion  strengths  of  2.4  and 
20.2  J/m2  respectively;  (c)  and  (d)  show  the  in-plane  tensile  stress  in  Y  direction 
siggYY  for  weak  and  strong  adhesion  strengths  of  2.4  and  20.2  J/m2  respectively. 
Larger  adhesion  strengths  lead  to  much  larger  tensile  stresses  in  the  sheet  ahead 
of  the  crack  bend.  Reduction  in  width  of  the  tear  leads  to  relieving  of  this  tensile 
stress,  and  thus  the  cracks  tapers  sharply  at  higher  strengths.  After  the  onset  of 
plastic  deformation,  load-transfer  at  the  polymer-graphene  interface  in  the 
composite  would  require  larger  contact  area,  and  a  weaker  interface  is  seen  to 
increase  the  size  of  the  tears,  thus  leading  to  larger  load  transfer. 
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M.J  Buehter,  MIT 


New  concept:  Biomimetic  silicon 


How  to  “transform”  a  brittle  material 
(here  as  model  selected  silicon)  into  a 
ductile  material. . . 


MwJL 


Diatoms:  Nature’s  flexible  armor 

Diatoms  are  unicellular  algae  with  silicified  cell  wall^ J  Buehler  MIT 
(frustules) 

-  Etymology:  dia  =  "through"  +  temnein  =  "to  cur.  Le.,  "cut  in  half rr 

-  One  of  the  most  common  types  of  phytoplankton 

Mechanical  properties  (using  FEM  /  experimental)1 

-  F.  kerguelensis  cells  of  50  mm: 


-  £, 


failure 


=  2.5%,  E  =  22.4GPa 


Raphe  alii 


FRU&TULE 

1.  Christian  E  Hamm  ef  at.,  Architecture 
and  material  properties  of  diatom  shells 
provide  effective  mechanical  protection. 

Atefune42i.  841^843  f20  February  2003} 

2:  Faldalore.  A.,  and  Bowler.  C.  REVEALING 
THE  MOLECULAR  SECRETS  OF  MARINE 
DIATOMS,  Annu.  Rev.  Plant  Biol.  2002.  53:109-30 


EPITHECA 


GIRDLE  BANDS 


Phaeoplast 


Golgi 


Overview  of  natural  biological  system,  which  shows  the  nanostructural  arrangement 
of  a  brittle  material  (silica). 

Can  we  interpret  this  design  and  use  it  for  the  design  of  new  composites? 
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Structural  hierarchies  in  two  different  silica-based  skeletal  assemblies  in  a  marine 
diatom  species  (a-d)  and  a  sea  sponge  (e-h).  The  diatom  in  (a)  Concinodicus  sp. 
has  a  silica-based  exoskeleton  (called  frustule)  made  up  of  porous  parts  arranged  in 
a  hierarchical  fashion  (a)  shows  the  whole  frustule  member  (external  surface  of  the 
diatom),  (b)  shows  areola  pores,  the  internal  surface  of  the  diatom(c)  shows  the  2nd 
central  porous  layer,  the  cribrum  and  (d)  shows  the  cribellum,  the  external  porous 
layer,  (e)  shows  the  external  cage  structure  of  the  silica-based  skeletal  system  of 
Euplectella  sp.  Scale  bar  5  mm  (f-g)  show  some  of  the  underlying  hierarchical 
structures  with  (f)  showing  fibre-composite  structure  in  a  constitutent  beam 
consisting  of  many  spicules,  Scale  bar  20  microns,  (g)  single  spicule  showing 
laminated  silica-protein  structure,  Scale  bar  5  microns  and  (h)  biosilica  constituent 
of  the  silica  layers  revealing  its  consolidated  nanoparticle  nature,  Scale  bar  500  nm. 
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Methodology 


M.J.  Buehler,  MIT 


•  Construct  bioinspired  design  of  silicon  nanostructure 

•  Load  in  tension  at  1 E-5  fs_1 

•  Temperature  300  K,  use:  ReaxFF-Lamrmps 


Brick  schematic 


silicon 


Periodic  boundary  condition  (non  fixed  ends) 


Equivalent  brick 


We  have  designed  a  first  model  system. 

The  following  studies  are  focused  on  varying  the  width  of  the  bricks  and 
testing  the  influence  on  the  mechanical  behavior. 
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Stress  fields  after  failure 


Results  show  that  the  larger  systems  feature  a  stress  concentration  at  corners; 
smaller  systems  deform  homogeneously. 
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Results  show  that  the  smaller  systems  behave  in  a  ductile  fashion,  and  feature 
large  deformation  with  >70%  strain. 
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Hierarchical  model  of  silica  material 


M.J.  Buehler,  MIT 


0  0,5  1  1,5  2  n  0  2  0  4  os  oo 


(a)  shows  the  geometry  of  the  nano-honeycomb  used  as  building  blocks  for  the 
composite  structures,  (b)  shows  a  section  of  the  triangular  mesh  mesoscale 
particle-spring  model  setup,  (c)  show  stress-strain  curves  obtained  from  atomistic 
simulations  of  a  nano-honeycomb  structure,  and  for  bulk  silica  with  a  crack  of  the 
same  size  as  the  pores  in  the  nano-honeycomb.  The  legend  shows  the 
classification  of  the  nano-honeycomb  structure,  which  is  shown  as  (t,  pi,  pw) 
parameters  for  the  structure,  values  given  in  A.  The  bulk  silica  structure  shows 
purely  brittle  fracture,  the  nano-honeycomb  structure  show  ductile  fracture,  (c) 
shows  the  behavior  of  the  mesoscale  triangular  mesh  lattice  fitted  to  this  constitutive 
behavior. 
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M 

J*$.  Hierarchical  material  design 


M.J.  Buehterr  MIT 


Develop  two 
design 
suggestions, 
where  the  same 
material  is  used 
in  the  entire 
structure 


However,  with 
different 
hierarchical 
distributions 


Randomly  distributed  particle  composite  structures  at  the  mesoscale.  Constituents 
are,  bulk  silica  (in  grey/light,  with  a  high  volume  fraction)  and  nano-honeycomb 
structures  (in  blue/dark,  small  volume  fraction).  Design  conditions  that  enhance 
toughness  of  bulk  silica  by  distribution  of  small  amount  of  nano-honeycomb  silica 
are  being  investigated  here.  The  structure  on  the  left  shows  bulk  silica  particles  as 
reinforcement,  whereas  the  right  structure  shows  nano-honeycomb  silica  particles 
as  the  reinforcing  phase.  Structures  of  both  types  are  studied  for  crack  propagation 
response. 
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(a)  i6o 
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Resulting  stress-strain  curves 


M.J  Buehtec  MIT 


strain 


strain 


Stress-strain  curves  for  (a)  bulk-silica  reinforcing  composite  structures,  with  and 
without  presence  of  a  pre-crack.  The  near-identical  response  shows  the  flaw- 
tolerance  behavior  for  these  structures  to  pre-cracks  of  the  given  size.  The 
structures  show  multiple  cracking  throughout  the  material,  and  this  is  reflected  in  the 
stress-strain  curve  as  a  gradual  loss  of  stiffness  of  the  material  as  the  number  and 
size  of  the  multiple  cracks  grow,  (b)  Stress-strain  curves  for  nano-honeycomb 
reinforcing  composite  structures,  with  and  without  presence  of  a  pre-crack.  The 
varying  fracture  strengths  clearly  show  an  effect  of  the  crack  size.  All  structures  fail 
by  the  growth  of  a  single  dominant  crack. 
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Deformation  &  failure  mechanisms 


0.48%  strain 


MJ  Buehlerr  MIT 


(b)| 

/,  \ 

/  1 

V 

*  1  f 

S 

0.64%  strain 


0.48%  strain  0,64%  strain 


Crack  pathways  for  composite  structures  with  nano-honeycomb  structure  as  the 
matrix  and  brittle  silica  as  the  reinforcing  particulate  phase.  The  volume  fraction  of 
silica  phase  is  76%.  (a)  and  (b)  show  fracture  progress  starting  from  a  material  with 
no  pre-crack;  (c)  and  (d)  show  fracture  progress  in  the  same  material  with  a  pre¬ 
crack  present.  In  cases  we  observe  that  the  pre-crack  propagates  for  a  small 
distance  but  does  not  propagate  through  the  sample,  and  other  smaller  cracks  are 
initiated  throughout  the  sample.  These  multiple  small  cracks  determine  the  stress- 
strain  response  of  the  structure.  The  structure  is  thus  flaw-tolerant  to  pre-cracks  of 
these  sizes,  and  the  fracture  stress  and  behavior  are  almost  independent  of  the  size 
of  the  pre-crack.  This  is  reflected  also  in  the  stress-strain  curve  for  the  stress-strain 
response  of  the  undetected,  and  cracked  structures  (shown  in  Figure  4a).  A  fracture 
toughness  cannot  be  measured  for  these  structures  and  crack  sizes. 
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Deformation  &  failure  mechanisms 


0.29%  strain  0.48%  strain 


Crack  pathways  for  composite  structures  with  brittle  silica  as  the  matrix  and  nano¬ 
honeycomb  structure  as  the  reinforcing  particulate  phase.  The  volume  fraction  of 
bulk  silica  phase  is  86%.  (a)  and  (b)  show  fracture  progress  starting  from  a  no  pre¬ 
crack  material;  (c)  and  (d)  show  fracture  progress  from  the  same  material  with  a 
pre-crack.  In  both  cases  we  observe  that  fracture  occurs  through  the  propagation  of 
a  dominant  crack,  (a)  and  (d)  show  the  un-cracked  and  pre-cracked  specimens  at 
the  same  load,  the  un-cracked  specimen  is  intact,  whereas  the  pre-crack  has 
started  propagating  in  the  other  specimen.  Since  the  fracture  strength  of  a  structure 
with  a  dominant  propagating  crack  is  pre-crack-size  dependent  (according  to 
fracture  mechanics),  the  stress-strain  curve  for  the  stress-strain  response  of  the 
undetected,  and  cracked  structures  are  markedly  different  (shown  in  Figure  4b). 
Fracture  toughness  can  be  measured  for  these  structures,  as  the  energy  required 
for  the  growth  of  the  pre-crack. 
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Different  composite  structures  with  brittle  bulk-silica  as  the  matrix  (volume  fraction 
86%)  and  ductile  nano-honeycomb  structures  as  the  reinforcing  particulate  phase 
showing  fracture  toughness  improvement  mechanisms.  The  particles  have  circular 
cross-section  and  are  randomly  distributed  and  five  different  random  structures  are 
shown  here.  A  single  pre-crack  is  introduced  and  then  subjected  to  mode  I  loading. 
Propagation  of  the  single  dominant  crack  is  seen  on  loading,  and  the  propagation 
path  is  marked  in  white.  All  structures  show  that  the  crack  path  is  not  straight,  but 
connects  reinforcing  particles  lying  close  to  the  original  crack  plane.  Crack 
deflection  and  bridging  by  reinforcing  particles  behind  the  crack  tip  are  the 
mechanisms  seen  to  increase  toughness  here. 
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J-integral  analysis  &  Ashby-plot 


M.J.  Buehler,  MIT 


(b) 


Crack  extension,  \a  (run) 


600% 

500% 

400% 

300% 

200% 

100% 


hierarchical  silica 
composite 


bulk  silica  i 


0% 


0%  20%  40%  60%  60%  100%  120% 


percent  change  in  modulus  (%) 


Left:  Calculation  of  the  J-integral  and  R-curves.  (a)  shows  the  J-integral  calculation 
for  a  stationary  crack  by  the  use  of  a  domain-integral  around  the  crack.  The  J- 
integral  provides  the  value  of  the  energy  release  rate  per  unit  advance  of  the  crack 
into  the  crack  tip,  or  the  resistance  to  crack  propagation.  The  red/dark  region  shown 
is  the  domain  of  integration  and  the  convergence  of  the  J-integral  is  tested  by  taking 
different  r,  and  r2  regions  for  the  same  crack  and  specimen  configuration. (b)  shows 
fracture  toughness  measures  as  a  function  of  crack  advance  (R-curve  behavior)  for 
all  the  structures  in  Figure  8.  The  toughness  of  bulk  silica  is  also  shown  as  a  dotted 
curve.  Right:  The  design  of  the  hierarchical  silica  composite  improves  the  toughness 
of  bulk  silica  significantly  (=5.5  times)  while  compromising  on  the  stiffness  only 
slightly  (=70%  of  bulk).  This  points  towards  the  use  of  hierarchies  along  with  a 
single  design  material  to  improve  undesirable  mechanical  properties  significantly 
(here  low  toughness)  while  not  compromising  on  the  desirable  ones  (here  high 
stiffness). 
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Q  Differential  multi-scale  modeling  of 
/SV  chemically  complex  materials  under  heavy 

# 


Today 


~ 100  lb. 
load 


deformation 

Optimization  of 
structure  at  multiple 
scales, 

simultaneously 


M.J.  Buehler,  MIT 
Objective  force 


<  30  lb. 
effective 
load 


How? 

Greater  material 
complexity,  merger 
70+  material  &  structure 


0  mph 


New  materials, 

reduced  weight 

at  higher  >  40  mph 

performance 


Adapted  from:  T.L.  Doligalski  (Engineering  Sciences  Directorate,  ARO) 


Our  research  suggests  new  opportunities  to  meet  challenges  of  a  lighter, 
faster  and  more  adaptable  Army  by  introducing  greater  material 
complexity.  Thereby,  the  use  if  hierarchies  and  novel  material  concepts  as 
studied  in  our  project  may  be  critical  to  achieve  the  target  properties. 

As  shown  earlier,  the  merger  of  structure  and  material  is  believed  to  be  a 
key  aspect  that  will  enable  us  to  meet  the  challenges  of  Defense 
applications. 
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Mir 


Massachusetts  Institute  of  Technology 


“Material  (pie)  in  the  sky” — what  if  we  could... 


■  ...  make  a  diversity  of  functional  materials  out  of  the  same 
universal  building  blocks? 

E.g.  use  simple  chemistry,  abundant  materials  (protein,  silica) 

■  ...make  materials  at  low  energy  cost,  with  low  weight  (density), 
and  yet  superior  properties  (all  of  the  above)? 

E.g.  avoid  high  temperature  processing 

■  ...make  adaptable  materials,  which  change  properties  based  on 
need,  e.g.  via  external  stimuli  (light,  electric  field,  pH,  water...)? 

E.g.  mechanical  response  changes  dramatically  based  on  threat 
type  (“modes”  of  operation) 


Have  in  abundance:  natural  materials 


water  &  rocks 


sand 


soy  beans  (protein) 


Need:  functional  materials,  devices,  armors 


"exaskeleton./ 
not  only 
provides  the 
protection  of 
traditional 
body  armor*  it 
enhances  the 
st  rength  of 
soldiers'  lops. 


RoboTroop 

Hie  US  Aimy^  vision 
for  2030.  Many  of 
these  technologies 
are  already  under 
developirtent. 


Unlocked  by  voice 
command  or 
electronic  tripper* 
the  weapon  would 
shoot  high  explo¬ 
sive  munitions  up 
to  IjOOG  meters, 
Another  feature; 

It  could  be  scaled 
from  non  lethal 
(’“stun")  to  lethal 
force. 


The  custonv 
fit  boot  will 
be  designed 
to  minimize 
effort  and 
Increase 
endurance. 


The  head gea r  has 

biometric  facial  recogni¬ 
tion  to  Identify  insur¬ 
gents*  while  targets  are 
1 1 1  umlnatcd  o  n  the  d  isplay. 


Advanced 
iianomaledal 
armor  protects 
against  blasts* 
burns  and 
rifle  rounds. 
The  entire 
ensemble  is 
embedded 
with  behav¬ 
ioral  and 
physiological 
sensors  that 
continuously 
monitor  the 
soldier's 
health. 


The  "data 
glove"  may 
housed 
to  operate 
robots  and 
unmanned 
drones. 
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SHIRK-  US  AflTfl 


Core  question:  How  to  turn  abundant  materials 

into  functional  materials? 


Less  energy,  lighter,  stronger,  adaptable, 
easily  available... 

Path  to  success  -  understand  biological 
materials 


Biological  highly  functional  materials  - 
diatoms,  sea  sponges,  organisms 
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Protein  materials  in  biology 


Eye’s  cornea 

(optical/mechanical, 
collagen  material) 


Skin  (protection, 
complex 
composite  of 
collagen, 
elastin) 


Nerves  (signaling) 


Cells  (complex 
material/system 
based  on  structural 
proteins 


Blood  vessels  (mass 
transport) 


Tendon 

(links  bone, 

Muscles,  “glue”,  e.g. 
Achilles  tendon,  collagen) 


Cartilage  (reduce 
friction  in  joints) 


Bone  (structural 
stability, 

collagen/minerals) 


Muscle  tissue 

(mechanical  force,  motor 
proteins) 


http://www.humanbodv3d.com/  and  http://publications.niqms.nih.gov/insidethecell/imaqes/ch1  cellscolor.jpg 


Genes  (DNA  structure)  define  protein 
structure 


ACGT  - > 

Four  letter 
code  “DNA” 

Transcription/ 

translation 


. .  -  Proline  -  Serine  - 
Proline  -  Alanine  -  . . 


Sequence  of  amino  acids 
“polypeptide”  (ID 
structure) 


Combination 
of  3  DNA 
letters  equals 
a  amino  acid 

E.g.:  Proline  - 

CCT,  CCC, 
CCA  CCG 


Mir 


Folding 
(3D  structure) 
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M.  Buehler  et  a/.,  Nature  Materials,  2009 


Protein  material  assembly 


Folded 

protein 


dimers 


intermediate  filaments 


axial  growth 


t 

• _ _ m 

■  *1. 


radial  growth 


1 


1 

* 


Self-assembly 

processes 

(e.g.  lamin  intermediate 
filaments) 


Experimental  analysis 

(Herrmann,  Aebi,  Kreplak 
et  al.) 


How  do  fiber  and  larger-scale  material  properties  (s.a. 
elasticity,  deformability,  failure,  fracture,  etc.)  depend  on  the 
structure  of  the  constituting  protein  structure? 
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M.  Buehler  et  al.,  Nature  Materials,  2009 


Tendon,  bone 
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Collagenous  tissues  -  an  issue  (tissue)  of  multiple 
scales 


Collagen  molecule 
length  s  300  nm 
diameters  1.6  nm 


macro 


Tissues  (tendon, 
bone) 

dimension  s  cm 


Collagen  fiber 
length  s  mm 
diameters  10  pm 


Collagen  fibril 
length  s  pm 
diameters  100  nm 


Amino  acids 

covalent  bond  si  A 

nano 


Bone 
Tendon 
Cartilage 
Eye’s  cornea 


M.J.  Buehler,  Proc  .  Nat'l  Academy  of  Sciences  USA,  2006;  R.  Ritchie,  M.  Buehler,  P.  Hansma,  Physics  Today,  2009 


Biological  materials  are  adaptable,  mutable... 


Biological  systems 
universally  feature  active 
components;  e.g.:  bone 
remodels  according  to 
changes  in  loading,  cells 
change  their  shape  under 
different  environments  or 
stress 

Mechanisms;  Sensing  of 
environmental  cues 
followed  by  gene 
regulation/activation,  which 
changes  the  hierarchical 
structure  of  materials 
(“mutability”) 


A 


New  blood  vessels  form  under  cyclic  strain 
to  repair  wounds,  provide  nutrients  where 
process  needed 


t 

Structure 


Property 

4 


1 


B 


/ 


Requirement 

Process  HN+* . 

Process  Hi 


Process  HI 

(transcription, 

translation) 


no  strain 


cyclic  strain 


Property  HN 
\  | 
Property///  / 


Structure  HO 


Structure  /// 


Structure  HN 


Materiomics 


Property  HO  . 


Gene  regulation/ 


activation 

Mechanotransduction  mechanisms  (e.g.  angiogenesis) 


'^Requirement  ftN 
Requirement  RJ 
Requirement  RO 


Sensing  (e  g.  charged 
requirements) 


Buehler  and  Yung,  Nature  Materials,  2009 


Y.  Yung,  M.  Buehler,  D.  Mooney,  et  al.,  PNAS,  2009 


Challenge:  Time-scales  on  the  order  of  minutes,  hours,  days 


Linking  chemistry,  structure  and  mechanics 


Molecular  dynamics: 

Fundamental  model  of  materials 
“bottom -up” 


rj 

^  system 


^  bond  ”1”  ^-  'finsrlp  +  j 


angle 


torsion 


“1“  ^  Coulomb  ^  vdW  + 


bond  stretching 


angle  bending 


bond  rotation 


covalent  bonds 


electrostatic 

interactions 


electrostatic  weak  bonds 
bonds 


vdW  interactions 


M.  Buehler,  Atomistic  Modeling  of  Materials  Failure,  Springer,  2008 
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Integration  with  experimental  techniques 


Time  scale 

A 


mm 


MS 


ns 


Nano- 

indenlation 


Micropipette 


Chemistry  - 

QM 


ps 


Meso- 

scale 

models 


Functionality 

MEMS  testing 


Continuum 

models 


Collagen  fibril 


Feed 

parameters 
Finer  trains 
coarser 


Oplical/magnetic  tweezers 


N  on- 
reactive 
MD 


Reactive 

MD 


(DFT) 


Tomography 


TEM  (e  g  cryo) 


x-ray 

diffraction 

ISIMR 


nm 


DNA 

polypeptides 


nanopartides 
( nanowires , 

CNTs)  ceils 


secondary  protein 
structures  |'SS,  AH, 
TC .  RC) 


tissues 

organs 

organisms 


m  Length  scale  Key  tool: 

Large-scale 

computing 


M.  Buehler  et  a!.,  Nature  Materials,  2009 


Pulling  on  a  single  collagen  molecule 


amino  acids 
-1  nm 


tropocollagen 
triple  helix 
-300  nm 


collagen 
fibrils 
~1  pm 

collagen  fiber 
fascicles... 
-10  pm 


M.J.  Buehler,  Proc .  Nat'l  Academy  of  Sciences  USA,  2006 
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Pulling  on  a  single  collagen  molecule 


Coarse-graining 

Mesoscale 


Mechanism: 

Entropic  elasticity 
(change  in  configurational 
entropy) 


Reduced  extension  (x/L) 


Advancement  in  experimental  equipment: 

Have  quantitatively  confirmed  predictions  from  our  simulations 

Sun  et  al.,  J.  Biomechanics,  2004;  Buehlerand  Wong,  Biophys.  J.,  2007 
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Hierarchical  features  of  soft  collagenous  tissue  I'lii 


300  nm,  collagen  molecule 


-> 


67  nm,  D  period 


amino  acids 

~1  nm 


tropocollagen 

-300  nm 


fibrils 

~1  (jm 


fibers,  fascicles 

-10  pm 


continuum 

»100  (jm 


Atomistic  & 

molecular 

modeling 


FEM 

(continuum) 

modeling 


1200 


X(Angstrom) 

70  MPa  (a)  I - 
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J 

B.  Moran,  M.  Buehler  et  a!.,  Annals  in  Biomedical  Engineering,  2009 


Materials  failure  in  biology  and  engineering 


Rapid  aging  disease  progeria:  mechanical  failure  of  cell  nucleus 


Patient 


Fracture  in  celi's  nucleus 
created  under  mechanical 
deformation 


mutation 


v>Vv< 

sAm  > 


/f 


Failure  of  protein  molecules 
(genetic  defect) 


Airplane  crash:  failure  of  engineering  structure 


microcrack 
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Rat  kangaroo  kidney  epithelial  -  immunofluorescently  labeled  for  keratin 
(intermediate  filament) 


Vimentin  intermediate  filaments  (IF)  Mlii 


Organize  the  internal  structure  of  cells 

Provide  strength  &  very  extensible 
(300%)  -  source  unknown, 
mechanisms  unknown,  . . . 


Cytoskeleton 


Full  length 
filaments 


Unit-length 

filaments 

(8  tetramers) 


Tetramers 

(2  dimers) 

Dimer 

(elementary 
building  blocks) 

20 


intermediate 
filaments 


>240nm 


10  nm 


-60  nm 


-45  nm 


Janmey  et  al. 


deformation 


nano 


Intermediate  filaments  -  widely  found  in  nature  I'lii 


Chromosomes 


hair,  hoof 


Nuclear 

fibroblast  cells 
(make  collagen) 


cell  nucleus 

Chromatin 


Intermediate  filaments  dimer  structure 


Crystallized:  Part  of  1 A  and  2B 
(structure  not  known  of  other  parts) 
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head 


j  tail 
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hydrophobic  stripe 


(f>  =-  58‘ 

if/  =-  4T 


Amino  acid 
sequence 

(genetics) 


Final 


Validation 


structure 


Initial 


structural 


model 


Energy 
minimization 
&  equilibration 

(repeat  until 
convergence) 
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x  m  :  ik1 


iX  . 


Z.  Qin,  L.  Kreplak,  M.  Buehler,  PLoS  ONE,  2009 


Structure  prediction  and  functional  properties  I'lii 


Genetics 

. .  .GGLGGQGAGA 
A A A A AG  GAG  QG 
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Molecular  structure  of  the  dimer  and  tetramer  INii 
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Z.  Qin,  L.  Kreplak,  M.  Buehler,  Nanotechnology,  2009;  PLoS  ONE,  2009 


Stress  (MPa) 


Stretching  of  intermediate  filament  protein 


Engineering  strain 


Z.  Qin,  L.  Kreplak,  M.  Buehler, 
PLoS  ONE,  2009;  PRL,  2010 
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Larger-scale  protein  network  properties 


Kreplak  et  al.,  2006 


How  does  behavior  of  protein  network  at  larger  scale  depend  on 

the  protein  constituent’s  structure  ? 


Lamin  intermediate  filament  network  (cell 
nucleus) 
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full-atomistic  model 


coarse-grained  model 


0(10.. 100  nm) 
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Ackbarow  et  a!.,  PNAS,  2007;  Ackbarow,  Buehler  et  al.,  PLoS  ONE,  2009 


Coarse-g raining  approach 
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Movie  of  deformation 


Intermediate 

filament 


force  8 


http://web.mit.edu/mbuehler/www/slides/HG- 

ASME/movie  filament.mpq 


Z.  Qin,  L.  Kreplak,  M.  Buehler,  PLoS  ONE,  2009; 
Z.  Qin  et  al.,  unpublished 


Deformation  behavior  of  alpha-helix  network 
under  strain 


strain 


oooo 


force 


Mir 
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Ackbarow,  Sen,  Thaulow, 
Buehler,  PLoS  ONE,  2009 


Self-protecting  mechanism 


B 
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Ackbarow,  Buehler  et  al.,  PLoS  ONE,  2009 


Hierarchical  levels  and  associated  mechanisms 
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H4  4- 

.  JL _ ] 
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L — h 

L  network  Structural  transformation  of  crack-like  defects 

to  mitigate  stress  concentrations 

H2 


square  lattice 
unit  cell 


alpha-helix 

filament 


Distance  between  filaments  facilitates  extreme 
strain  gradients  at  low  energy  cost 

Parallel  AHs  &  long  linear  array  of  turns 
provides  structural  basis  for  large  extensibility 
via  repeated  rupture  of  turns;  A-B-transition 


Clusters  of  3-4  H-bonds  provide  optimal 
resistance  against  mechanical  failure  (3-4  H- 
bonds  break  concurrently) 


H-bond 


H-bonds  form  at  moderate  temperatures:  Self- 
assembly  of  alpha-helices/self-healing 


Ackbarow  et  a!.,  PNAS,  2007;  Ackbarow,  Sen,  Thaulow,  Buehler,  PLoS  ONE,  2009 
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I  will  highlight  of  few  of  the  disciplines  in  which 
we  have  elected  to  increase  our  efforts.  In 
FY1 1 ,  the  Army  will  increase  its  focus  on 
developing  the  scientific  foundations  leading  to 
the  discovery  of  novel  materials  with 
extraordinary  performance  characteristics  of 
particular  interest  to  the  Army,  such  as  ballistic 
protection.  The  robust  research  approach  will 

emphasize  multi-scale  modeling  and 

simulation  for  performance  prediction  and 
design  of  materials  under  extreme  conditions 
(high  strain  rate,  temperature,  pressure, 
etc.).  This  modeling  effort,  supported  by  in- 
house  and  extramural  experimental  studies, 
will  study  a  range  of  material  classes  (such  as 
structural,  electronic,  energetic)  and  coupled 
disciplines  (e.g.,  mechanics  and 
electromagnetics)  enabling  two-way 
information  transfer  across  the  length  scales 
from  the  molecular  level  up  to  the  material 
system  level. 


Courtesy  J.  Andzelm  (ARL) 
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Mechanical  response  of  spider  silk 
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Z.  Shao  and  F.  Vollrath,  Nature,  2002 


Spider  silk  with  diameter  of  O(inch)  can  stop 
large  aircraft  in  flight 
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Keten  and  Buehler,  Nano  Letters,  2008,  Buehler  et  a!.,  Prog.  Mat.  Science,  2008 


Spider  silk’s  hierarchical  structure 
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core 
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protein  motif 
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H -bonded 
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Scaling  of  strength  of  H-bond  assemblies 


Characteristic  maximum  number  of  H- 
bonds  work  cooperatively:  3. .4  H-bonds 

___ _ _  ,  ,  — ► 

1 1 1  i  1 1 1  i  I  i  n  1 1  i  5 1  n  r 


Many  small  segments  with  3-4  HBs  each 
enhance  the  strength  multiple  times 


Intrinsic  energetics, 
nanomechanics,  and  structural 
confinement  control  strength 
properties:  Reach  strength  at 
characteristic  dimension  of  H- 
bond  clusters 

Extension  (nm) 

Keten  and  Buehler,  Nano  Letters,  2008,  Buehler  et  a!.,  Prog.  Mat.  Science,  2008 


Size  effects  of  strength  of  beta-sheet  proteins 
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Keten  and  Buehler,  Physical  Review  E,  2008 


Validation  against  experiment 


Simple  model  explains  many  experiments, 

has  implications  for  a  wide  range  of  polymeric  materials 
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Theory  predicts  near-equilibrium  strength  of  beta-proteins  accurately  to  be 
around  50-300  pN. 


Keten  and  Buehler,  PRL,  2008,  Ackbarow,  Keten  et  al.,  J.  of  Phys.  Cond.  Matt.,  2009 


characteristic  dimension 
(#  H-bonds) 


Cluster  size  in  natural  protein  materials 


Strength  limit 
(prediction) 


Alpha- helices 


Beta-sheets 


Beta- helix 


S.  Keten  and  M.J.  Buehler,  Nano  Letters,  PRL,  2008;  M.  Buehler  and  S.  Keten,  Rev.  Mod.  Phys.  2010 


Properties  of  H-bonds 


H-Bond  energy:  2-10  kcal/mol 

Thermal  energy  scale 
kBT~  0.6  kcal/mol  (room  temp) 


controls  unique  properties  of  water 

Water:  liquid  @  300K  ia/uy9 
Protein:  solid  @  300K  VVn  Y  ' 


individual  H-bonds  in  water 
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Hierarchical  assembly  -  from  nano  to  macro 
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filament 


alpha-helix 

turn 


Protein  filaments  and  protein  materials 
constitute  hierarchical  assemblies,  from 
nano  to  macro 

H-bonds  provide  a  “universal  glue”  found 
proteins  (e.g.  alpha-helices,  beta-sheets, 
etc.)  -  Nature’s  cement 

Fundamental  question — mechanical 
properties  of  H-bonds 


42  HO  •  ■  ■  OO  H-bond 


Universal  building  blocks,  diverse  structures  (and 
thus,  functions) 
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M.  Buehler  et  a!.,  Nature  Materials,  2009 
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Mechanical  response  of  spider  silk 
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Z.  Shao  and  F.  Vollrath,  Nature,  2002 


Spider  silk’s  hierarchical  structure 
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H-bond  (chemical 
structure) 

nano 


S.  Keten,  M.  Buehler, 
APL,  2010 


Structure  prediction  and  functional  properties 


Genetics 
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Replica  Exchange  approach  (parallelized) 


(temperature) 


LU 


E 


Simulate  copies  of  same  system  with  different  temperatures 
Exchange  configurations  between  them 

Evaluate  most  stable  configurations  obtained  at  low  temperature  = 

solution  “ensemble”  of  structures 
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S.  Keten,  M.  Buehler,  APL,  2010 


Search  for  most  likely  conformations 


Key  finding:  Secondary  structure  is  dominated  by  beta-sheets  (mostly 
seen  in  poly-Ala  regions),  as  well  as  turn  type  structures.  This  is  in  line 
with  what  has  been  observed  in  spider  silk.  (%  below  show  relative 
amount  of  respective  secondary  structure) 


HELIX:  0% 

BS  :  49.78% 
COIL  :  19.56% 
TURN  :  30.67% 


HELIX:  0.44% 
BS  :  40.44% 
COIL  :  15.67% 
TURN  :  43.44% 


S.  Keten,  M.  Buehler, 
APL,  2010 


Probability 


Experimental  validation 
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0.08 
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Replica  #42 


CA(i+3)  distances 
along  chain  =  9  A 


4  6  8  10 

Distance  CA(i)-CA(i+3)  (A) 


experiment  simulation 


•  o  « n 


alanine  dihedral  angle  distribution 


-  poly-Ala  regions  form  beta-sheet  regions  with  stacking  to  accommodate 
methyl  side-chains.  These  interlocked  regions  are  densely  packed  and 
orderly. 

-  Glycine  residues  predominantly  lie  in  the  31  helix  structural 
conformation,  with  occasional  type  II  beta-turn  structures.  31  helix 
structure  is  supported  by  CA(i)-CA(i+3)  distances,  being  most  probable 
around  9  A,  as  well  as  existence  of  interchain  H-bonds.  Packing  of  chains 
is  disorderly. 
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S.  Keten,  M.  Buehler,  APL,  2010 


Detailed  view  of  structure 
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S.  Keten,  M.  Buehler,  APL,  2010 


Mechanical  analysis  of  spider  silk 


Mir 


■  Apply  force  to  half  of  the  strands  in  each  direction  (7  total).  Strands 
are  randomly  selected  to  mimic  loading  conditions  in  spider  silk. 

■  A  constant  force  is  applied  to  each  strand  and  systems  is 
equilibrated.  Loading  rate  (pN/time)  is  kept  constant  for  each  addition 
of  force  [here  0.2  pN/ps] 

■  Extension  is  taken  from  the  center  of  mass  of  pulled  atoms, 
converted  to  strain  based  on  initial  distance. 

■  Cross  sectional  area  taken  as  4  nm2  for  all  systems 


S.  Keten,  M.  Buehler,  APL,  2010 


Force-deformation  behavior 


End-to-End  Length  (nm) 


S.  Keten,  M.  Buehler,  J.  Roy.  Soc.  Interface,  2010 


Deformation  analysis  of  spider  silk 


S.  Keten,  M.  Buehler, 
J.  Roy.  Soc.  Interface, 

2010 
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Keten,  Xu,  Ihle,  Buehler,  Nature  Materials,  2010 
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Size  effects  in  silk  nanocrystals 


pull-out  setup  (small  system) 


stick-slip  mechanism 


pull-out  setup  (large  system) 


crack  formation  &  propagation 


Keten,  Xu,  Ihle,  Buehler,  Nature  Materials,  2010 


Size  effects  in  silk  nanocrystals 
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Keten,  Xu,  Ihle,  Buehler,  Nature  Materials,  2010 


Size  effects  in  silk  nanocrystals 


Interplay  of  scaling  laws  in  spider  silk 
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Natural  silk:  L*= 2-7  nm 
and  h*  =2-3  nm 


Keten,  Xu,  Ihle,  Buehler,  Nature  Materials,  2010 


Larger-scale  model  of  silk  (1000  nm) 


b 

1500 


Displacement  (AI*)  [A] 


F 


F 


- - > 

TTTTTT 


How  do  nanoscale 
structures  &  properties 
define  larger-scale 
mechanical  properties? 


Keten,  Nova,  Buehler  eta!.,  Nano  Letters,  2010 
Nova  et  al.,  in  preparation 


Stress  [MPa] 


Micro/mesoscale  mechanics  of  spider  silk 


Simulation 


Molecular  model  explains  experimental  results:  Size  effect  a  direct  result  of 
structural  changes  with  profound  impact  on  macroscale  mechanical  behavior 
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Keten,  Nova,  Buehler  et  al.,  Nano  Letters,  2010;  Nova  et  al.,  in  preparation 


Merger  of  structure  and  material 


Macroscale  structural  engineering 
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Ultra  scale  structu  ral  engineer!  ng 

Buehler  et  a!.,  Materials  Today ,  2007 


Universality  and  diversity  (universality-diversity 
paradigm) 
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(diversity)  by 
changing 
structural 
arrangements  of 
few  (universal) 
constituents 
rather  than 
inventing  new 
building  blocks 
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multifunctionality 


M.  Buehler,  Nature  Nanotechnology,  2010 


Concept  -  mutable  materials 


Concept:  Change  structure  of  material  at  defined  scales 
Thereby:  Induce  change  in  mechanical,  optical,  thermal  etc. 

properties 
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M.  Buehler,  Nature  Nanotechnology,  2010 


Practical  impact:  Transfer  to  synthetic  fiber 
design  based  on  nanomaterials 
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Wood-inspired  CNT  bundles 
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Bioinspired  thermal  materials 


heat  source 


Funded  by 
DARPA-MTO 


heat  source^ 


Inspired  by  biological  structures 

(e.g.  cell,  left),  controlled  assembly 
of  nanostructures  into  functional 
materials  (e.g.  by  using 
hierarchies),  some  of  the  challenges 
can  be  addressed 

Xu  and  Buehler,  Nano  Letters,  2009 
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H’osrtiori  x 


diatoms  (silicified  algae) 


Hierarchical  structure  of  sea  sponges,  diatoms  I'lii 


Hierarchical 
structure  of 
diatoms,  showing 
their  porous  silica 
structure 


Garcia  and  Buehler,  COMMAT,  2010  v — v— J 
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Silica  and  silicon  is  abundant,  but  brittle 
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most  abundant 
material  in  earth’s 
crust 
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Sen,  Buehler,  eta!.,  PRL,  2010 


Diatoms  -  Nature’s  flexible  armor 
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Garcia  and  Buehler,  COMMAT,  2010 
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Mechanism  of  deformation 
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Garcia  and  Buehler,  COMMAT,  2010 


An  opportunity  for  the  future 


I'lii 
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Tank  &  armor  out  of  glass  &  soy 
bean? 

Abundant,  cheap  (low  energy), 
stores  C02  (diatoms)... 


Summary,  conclusion  &  discussion 


■  Big  computers  allow  for  atomistic  simulation:  New  paradigm  of 
studying  protein  materials  -  “bottom-up”  approach 

■  Protein  materials  must  be  understood  at  multiple  scales  -  from  the 
level  of  individual  proteins,  to  fibrils,  fibers,  tissues,  etc. 

■  Mechanical  properties  matter  in  biology  -  understanding  failure  is  the 
key  to  success: 

Source  of  strength,  robustness,  adaptability  (remodeling) 

Mechanisms  in  diseases  and  injuries  (genetic  diseases,  cancer, 

etc.) 

Transfer  to  bioinspired  designs  (engineering) 
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